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Abstract 

For two-tier networks consisting of macrocells and femtocells, the channel access mechanism can 
be configured to be open access, closed access, or hybrid access. Hybrid access arises as a compromise 
between open and closed access mechanisms, in which a fraction of available spectrum resource is shared 
to nonsubscribers while the remaining reserved for subscribers. This paper focuses on a hybrid access 
mechanism for multi-channel femtocells which employ orthogonal spectrum access schemes. Consid- 
ering a randomized channel assignment strategy, we analyze the performance in the downlink. Using 
stochastic geometry as technical tools, we derive the distributions of signal-to-interference-plus-noise 
ratios, and mean achievable rates, of both nonsubscribers and subscribers. The established expressions 
are amenable to numerical evaluation, and shed key insights into the performance tradeoff between 
subscribers and nonsubscribers. The analytical results are corroborated by numerical simulations. 
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Channel management, femtocell, hybrid access, spatial Poisson process, two-scale approximation, 
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I. Introduction 

In current cellular network services, about 50% of phone calls and 70% of data services 
take place indoors [lj. For such indoor use cases, network coverage is a critical issue. One 
way to improve the indoor performance is to deploy the so-called femtocell access points 
(FAPs), also known as home base stations, besides macrocell base stations (MBSs). Femtocells 
are small cellular base stations, typically designed for use in home or small business EOS. 
The use of femtocells not only benefits the users, but also the operators. As the distance 
between the transmitter and receiver is reduced, users will enjoy high quality links and power 
savings. Furthermore, the reduced transmission range also creates more spatial reuse and reduces 
electromagnetic interference. 

Among the many challenges faced by femtocells, and more generally, two-tier networks, is 
the issue of interference; see Figure [Q The two-tier interference problem differs from that in 
traditional single-tier networks in several important aspects: First, due to the limitations of access 
mechanism, a user equipment (UE) may not be able to connect to the access point which offers 
the best service. Second, since femtocells connect to operator's core network via subscribers' 
private ISP, coordination between macrocells and femtocells and among femtocells is limited. 
Finally, compared to planned macrocell deployments, femtocells are usually deployed in an ad 
hoc manner, and the thus randomly placed femtocells make it difficult to manage the interference. 
In two-tier networks comprising of macrocells and femtocells, interference can be categorized 
into two types: 

• Cross-tier: This refers to the interference from one tier to the other; that is, femtocell 
interference at a macrocell and macrocell interference at a femtocell. 

• Co-tier: This refers to the interference within a tier; that is, femtocell interference at a 
nearby femtocell and macrocell interference at a nearby macrocell. 

In this paper, we consider two-tier networks that are based on multicarrier techniques, for 
example those deploying LTE or WiMAX standards, which use orthogonal frequency-division 
multiple access (OFDMA) techniques. In multicarrier systems, the available spectrum is divided 
into orthogonal subcarriers, which are then grouped into multiple subchannels, assigned to dif- 
ferent users. Due to the flexibility in channel assignment, both cross-tier and co-tier interferences 
may be alleviated. 
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The access mechanism of femtocells (see, e.g., flU) is a key factor that affects the perfor- 
mance of two-tier networks, and generally can be classified as follows, where we call the UEs 
that register to a femtocell as subscribers, and those that do not register to any femtocell as 
nonsubscribers. 

• Closed access: An FAP only allows its subscribers to connect. 

• Open access: An FAP allows all its covered UEs, no matter registered or not, to connect. 

• Hybrid access: An FAP allows its covered nonsubscribers to connect via a subset of its 
available subchannels, and reserves the remaining subchannels for its subscribers. 

Hybrid access flU is an intermediate approach, in which a fraction of the total available resource 
is allocated to nonsubscribers. By doing so, nonsubscribers near an FAP may handover into the 
femtocell to avoid high interference; meanwhile, with a certain amount of resource reserved 
for subscribers, the performance of subscribers may be well assured even in the presence of 
nonsubscribers. 

In hybrid access, a central issue is how to allocate the resource between subscribers and 
nonsubscribers. Previous studies (61 indicate that hybrid access improves the network perfor- 
mance at the cost of reduced performance for subscribers, therefore suggesting a tradeoff between 
the performance of nonsubscribers and subscribers. In this paper, we consider a hybrid access 
mechanism that uses a randomized channel assignment strategy, and analyze the performance in 
the downlink of both macrocells and femtocells. We employ stochastic geometry to characterize 
the spatial distributions of users as well as access points; see, e.g., jH and references therein for 
its recent applications in wireless networks. Accordingly, we derive the key performance indi- 
cators including blocking probabilities, mean achievable rates, and distributions of the signal-to- 
interference-plus-noise ratios (SINRs), of both nonsubscribers and subscribers. In our study, we 
establish general integral expressions for the performance indicators, and closed form expressions 
under specific model parameters. With the obtained results, we reveal how the performance of 
subscribers and nonsubscribers trades off each other. 

An overview of related works is as follows. In [9|, the authors proposed to study key per- 
formance indicators for cellular networks, such as coverage probabilities and mean achievable 
rates, through models and tools in stochastic geometry. Our work explores that approach in hybrid 
access femtocell networks, focusing on multi-channel models and the downlink performance. In 
ifTOl . the considered scheme divides the spectrum resource into two orthogonal parts which are 
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assigned to macrocells and femtocells, respectively, with femtocells being closed access. In iTTlTl . 
the authors considered two-tier femtocell networks using time-hopped CDMA, examining the 
uplink outage probability and the interference avoidance capability. In Ifl2l and |fT3l , the authors 
studied the performance of various femtocell access mechanisms, under substantially different 
system models from ours. Our work examines the interaction among multiple MBSs and FAPs, 
and focuses on scenarios in which macrocells and femtocells share the same frequency band, 
following a hybrid access mechanism, as will be elaborated in Section UH 

The remaining part of the paper is organized as follows. Section [XT] describes the two-tier 
network model, the channel assignment strategy, and the hybrid access mechanism. Based on 
a two-scale approximation for the spatial distributions of femtocells, Section Hn] analyzes the 
statistical behavior of UEs, deriving the distributions of the number of UEs connecting to either 
an FAP or an MBS, as well as the probabilities of a subchannel being used by either an FAP 
or an MBS. Built upon those statistics, Section [IV] establishes expressions for the distributions 
of the SINRs, and the mean achievable rates. Section [V] illustrates the aforementioned analysis 
by numerical results, which are also corroborated by simulations. Finally, Section IVT1 concludes 
the paper. 

II. Network Model 

A. Hybrid Access Femtocells 

In the two-tier network, we consider two types of access points, MBSs and FAPs. The MBSs 
constitute the macrocell tier, and they induce a Voronoi tessellation of the plane (see Figure [2]). 
When a UE attempts to access the macrocell network, it chooses to connect to the MBS in the 
Voronoi cell which it is situated in. An FAP provides network access to UEs in its vicinity, and 
we assume that all FAPs have a covering radius of Rf. Within the covered circular area of each 
FAP are two types of UEs, called subscribers and inside nonsubscribers. Inside nonsubscribers 
are those UEs who gather around an FAP without subscribing to its service; for example, transient 
customers in a shop or restaurant. Besides those two types of UEs, we also consider a third type 
of UEs, outside nonsubscribers, who are uniformly scattered over the whole plane, corresponding 
to those regular macrocell network users. 

The available spectrum is evenly divided into M subchannels, which are to be shared by both 
the macrocell tier and the femtocell tier. Each FAP is configured to allocate a fixed number, M s , 
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of subchannels for its covered inside nonsubscribers. These M s subchannels are called shared 
subchannels, and the remaining M r = M — M s subchannels are called reserved subchannels 
as they are reserved for the subscribers. In the considered hybrid access mechanism, each FAP 
selects its shared subchannels randomly, and independently of other FAPs. We assume that 
each UE, whether subscriber or nonsubscriber, needs one subchannel for accessing the two-tier 
network. When a UE accessing an MBS or an FAP, the serving subchannel is selected randomly 
(see Figure [3]). 

The hybrid access mechanism operates as follows. 

• A subscriber accesses to one of the M r reserved subchannels of its corresponding FAP. 
When there are more than M r subscribers in an FAP, they are served by time-sharing with 
equal time proportion. 

• An inside nonsubscriber accesses to one of the M s shared subchannels of its covering 
FAP. When there are more than M s inside nonsubscribers in an FAP, they are served by 
time-sharing with equal time proportion. 

• An outside nonsubscriber accesses the MBS in the Voronoi cell which the outside nonsub- 
scriber is situated in. When there are more than M outside nonsubscribers in the Voronoi 
cell, they are served by time-sharing with equal proportion. 

B. Mathematical Model and Two-scale Approximation 

To formulate the aforementioned hybrid access scenario mathematically, we model the spatial 
distributions of the nodes using spatial Poisson processes as follows. The MBSs constitute a 
homogeneous Poisson point process (PPP) $ m of intensity A m on the plane; the FAPs constitute 
another homogeneous PPP $j of intensity A/. In the circular covered area of radius Rt of each 
FAP, the subscribers are distributed according to a homogeneous PPP of intensity A s , and the 
inside nonsubscribers are distributed according to another homogeneous PPP of intensity A in . 
Outside nonsubscribers constitute on the whole plane a homogeneous PPP of intensity A out . All 
the PPPs are mutually independent. 

In this paper, we focus on the downlink performance. The transmit power is set to a constant 
value P m for an MBS, and Pf for an FAP. For the sake of convenience, we adopt a standard path 
loss propagation model with path loss exponent a > 2. Regarding fading, we assume that the 
link between the serving access point (either an MBS or an FAP) and the served UE experiences 
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Rayleigh fading with parameter fi. The received signal power of a UE at a distance r from its 
serving access point therefore is P m hr~ a (MBS) or Pfhr~ a (FAP) where h ~ Exp(/i). The 
fading of interference links may follow an arbitrary probability distribution, and is denoted by 
g. Furthermore, considering the typical scenario of indoor femtocell deployment, we introduce a 
wall isolation at the boundary of each FAP coverage area, which incurs a wall penetration loss 
factor W < 1. For all receivers, the noise power is a 2 . 

The different PPPs corresponding to different entities in the network interact in a complicated 
way, thus making a rigorous statistical analysis extremely difficult. For example, an inside 
nonsubscriber may be covered by more than one FAPs, thus leading to the delicate issue of 
FAP selection, and furthermore rendering the subchannel usage distributions among FAPs and 
MBSs intrinsically correlated. To overcome the technical difficulties due to spatial interactions, in 
the subsequent analysis we propose a two-scale approximation for the network model, motivated 
by the fact that the covered area of an FAP is significantly smaller than that of an MBS. The 
two-scale approximation consists of two views, the macro-scale view and the micro-scale view. 
The macro-scale view concerns an observer outside the coverage area of an FAP, and in that 
view the whole coverage area of the FAP shrinks to a single point, marked by the numbers 
of subscribers and inside nonsubscribers therein. The micro-scale view concerns an observer 
inside the coverage area of an FAP, and in that view the coverage area is still circular with 
radius Rf in which the subscribers and inside nonsubscribers are spatially distributed. By such 
a two-scale approximation, an inside nonsubscriber can only be covered by a unique FAP, and 
the coverage area of an FAP can only be within a unique Voronoi cell of an MBS. These 
consequences substantially simplify the performance analysis. In Section |V] we validate the 
two-scale approximation method through comparing analytical results and simulation results, for 
network parameters of practical interest. 

III. Statistics of UEs and Subchannels 

In this section, we characterize the distributions of UEs connecting to different types of access 
points, and the distributions of used subchannels in MBSs and FAPs. The analysis is based on a 
snapshot of the network model, and the obtained results will then be applied for characterizing 
the distributions of SINRs and achievable rates in Section HVl 
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A. Distributions of UEs 

Let U s be the number of subscribers accessing a given FAP and from our model we have 
U s ~ Poisson(A s 7r.Rj). Similarly, let U in be the number of inside nonsubscribers accessing a 
given FAP, and we have U m ~ Poisson(Ai n 7ri?j). 

The number of outside nonsubscribers who access a given MBS, denoted by U ont , is charac- 
terized as follows. We note that the macrocell coverage area is a Voronoi cell, and denote by S 
the area of the Voronoi cell. There is no known closed form expression of the probability density 
function (pdf) of S, whereas a simple approximation 031 has proven sufficiently accurate for 
practical purposes. Considering scaling, the approximate pdf of the size of a macrocell coverage 
area is given by 

f(S) = ^/X(SA m )i exp(-^A m )A m . (1) 

15 V ^ Z 

Conditioned upon S, the number of outside nonsubscribers is a Poisson random variable with 
mean \ out S. The probability generating function (pgf) of the unconditioned U out is thus given 

by 

G(z) = J™ exp (\ on% {z - l)S)f(S)dS. (2) 

Plugging in the approximate pdf of S and simplifying the integral, we get 

343 /7/7 A out y\ 

The distribution of U ont is therefore given by the derivatives of G(z), 

PRut = t} = —r 1 , ^ = o, l, ... . (4) 

B. Distributions of Subchannel Usage 

Since the subchannels are uniformly and independently selected by each FAP, it suffices to 
analyze an arbitrary one of them. Let us examine the probability that a given subchannel is used 
by an MBS or an FAP. In order to obtain the probability, we evaluate the average number of 
subchannels used by an MBS or an FAP, and then normalize the average number by the total 
number of subchannels, M. 
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The probability that a subchannel is used by an FAP is 

^ oo oo 

Pbusy,/ = ^(Yl mil1 ^' M rM U * =0 + X] mil1 ^' M sMUin = j} 
i=0 j=0 

.M r M OO 

= = o + a - if) E p {^ = o + 



«=0 i=M r +l 
-.M s oo 

j=0 j=M s + l 

For a Poisson random variable iV ~ Poisson(A), its cumulative distribution function (cdf) is 
given by 

± m = (i = ±^ = *±hx>, 

i=0 i=0 

where T(s,x) = f°° i s_1 e~*dt is the incomplete gamma function. Using © to simplify -Ptmsy,/* 
we get 



Pbus y ,/ = 1 + (A s 7r J Rjr(M r , A^Pj) - T(M r + 1, A.ttPJ)) + 

+ ^T4T(A i „7r J Rjr(M s ,A in 7ri?J) -r(M s + l,A in 7ri?J)). (7) 



M M s 

The probability that a subchannel is used by an MBS is 

j oo 

n usy , m = — E min {^^}p{^out = ^}, (8) 

i=0 

where P{C/ out = z} is given by ©. 

The spatial process of FAPs that use a given subchannel is the independent thinning of the 
original PPP of FAPs $j by the probability -Pbusy,/* denoted by with intensity A^ = XfPbusyj- 
The term "independent thinning" means that $^ can be viewed as obtained from by inde- 
pendently removing points with probability 1 — Pbusy,/- Similarly, the spatial process of MBSs 
that use a given subchannel is the independent thinning of the original PPP of MBSs $ m by the 
probability .Pbusy,™, denoted by <3>' m with intensity X' m = A m Pbu Sy ,m- These two independently 
thinned PPPs will prove useful in the subsequent analysis. 

IV. Analysis of SINR and Mean Achievable Rate 

In this section, we derive the distributions of SINRs and the mean achievable rates of all the 
three types of UEs. For each type of UEs, we begin with general settings, and then simplify the 
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general results under several specific parameters to gain insights. The mean achievable rates are 
the averaged instantaneous achievable rates over both channel fading and spatial distributions of 
UEs and access points. 

A. Macrocell UEs 

1) General Case: For an active UE served by an MBS, it must be occupying one subchannel 
of the MBS. The following theorem gives the cdf of SINR and the mean achievable rate of each 
active macrocell UE in general case. 

Theorem 1: The cdf of the SINR of a macrocell UE, denoted by Z m (T) = P{SINR < T}, 
is given by 

Z m (T) = 1 - nxj exp (nv\' m (l - /3(T, a) - -^—) - (9) 

Jo ' busy,m m ' 

The mean achievable rate of a macrocell UE is given by 

exp (nv\' m (l - /3(e* - l,a) - ) - ^ ± } -)dtdv. (10) 

^ busy,m * ra ' 

In © and fl©, (3{T, a) is given by 

ma) = ^ E Li(r ( - 2 -,,T g) - (i + 

The proof of Theorem [T] is in Appendix [A] 

In Theorem [H Z m (T) in © gives the probability that the SINR is below a given target level 
T, and r m in (TTOl) gives the mean achievable rate of a macrocell UE. The integrals in © and 
(flOl) can be evaluated by numerical methods, and furthermore, they can be simplified to concise 
forms in a number of special cases. The cases of interest in this paper are the combinations of 
(i) interference experiencing Rayleigh fading, (ii) the path loss exponent being a = 4, and (iii) 
no noise or the noise is weak, a 2 — > 0. The following discussions are categorized into two parts 
according to the interference fading. 

2) Special Case When Interference Experiences General Fading: 

a) General Fading, Noise, a = 4: We can simplify © via invoking the identity 
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where Q(x) = h= J r ^° exp(— y 2 /2)dy is the standard Gaussian tail probability. Letting a(T) 
7rA' m (/3(T,4) + X m /X' m - 1), we have 



ZJV = 1 - ,A ml /^exp (^) Q ( a(r) ^). (13) 

For a given distribution of interference fading, j3(T, 4) can be evaluated by numeric methods, 
and subsequently the SINR distribution Z m (T) can be obtained conveniently. 

b) General Fading, Small Noise, a > 2: Consider the interference-limited regime in which 
the noise power is approaching zero. Using the expansion e~ x = 1 — x + o(x) as x — > 0, the 
SINR distribution behaves like 



2\ 



m( } = " 1 + (B(T a) - 1) P h + 7 / xx«/2+i + °( a 

(14) 

When a 2 = 0, the SINR distribution is further simplified into 

Zm ^ = 1 - ,+(r(tJ up • (15) 

1 + {p[l,a) - l)-Pb US y, m 

3) Special Case When Interference Experiences Rayleigh Fading: Here we consider the case 
where the interference experiences Rayleigh fading with mean i.e. g ~ Exp(/i). In this case, 
the results are simplified as follows. 

Corollary 1: When the interference follows Rayleigh fading, the cdf of the SINR is 

Z m (T) = 1 - n\ m J exp ( - irv\ m - irv\' m ip(T, a, T~ 2/a ) - 

fPfWT \ ^V\. 
_7 ™ A M p — p ) dv - ( 16 ) 

The mean achievable rate is 

= ^X m J J exp ( - 7ruA m - irvX^ip^e* - I, a, (e* - l) _2/o ) - 

—7Tv\f(p ( ^ i,a,0j Mdudt, (17) 



7~n 



where 



The results in Corollary Q] are further simplified in the following special cases. 



(18) 
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a) Rayleigh Fading, Noise, a = 4: Using (fT2l) . we get 



Zm{T) = 1 - .A ra ,/^exp (^) Q ( K (T)^). (19) 
The mean achievable rate of a macrocell UE is simplified into 



'/^^^ ex KC^^)^(^'-Vw^y )d ' <20) 

where 

k(T) = Ti\ m + 7rA' m VT(7r/2 - arctan VT) + p ■ (21) 

^ V -*ra 

Z?) Rayleigh Fading, No Noise, a > 2: When a 2 = 0, the SINR distribution is simplified 

into 

Z m (T) = 1 -, . (22) 



The mean achievable rate of a macrocell UE is simplified into 

1 



1 + ^(e* - 1, a, (e* - !)-*/«) + ^>(™?^, a, 0) 



dt. (23) 



Specifically, when a = 4 we obtain 



z ™ ( T ) = 1 7 v TwWS • (24) 

1 + y/T{* - arctan VT + z^*£l)P busy>m 



The mean achievable rate of a macrocell UE is simplified into 

2 



tan?/ + (tt/2 - y)P b u Sy ,m + f jjr- 



:dy. (25) 



From a practical perspective, it is desirable to shape Z m (T) to make it small for small 
values of T. From ([24]) . we see that there are two approaches to shape Z m (T). First, Z m (T) 
decreases as Pbus y ,m, the probability that a subchannel is used by an MBS, decreases. This may 
be interpreted as an effect of frequency reuse. Second, Z m (T) decreases as the whole term, 
y-yW^^ decreases, which corresponds to a number of network parameters, representing the 
effect due to the deployment of the femtocell tier. 
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B. Femtocell UEs 

1 ) General Case: A UE served by an FAP occupies one subchannel of the FAR The following 
theorem gives the cdf of SINR and the mean achievable rate of each active femtocell UE in 
general case. 

Theorem 2: The cdf of the SINR of a femtocell UE in general case is given by 

Z f( T ) = 1 -Jp f l ex P ( - P{oc)T 2/a v - * ^ )dv. (26) 
The mean achievable rate of a femtocell UE is given by 

V - 4 f ? f - P ( - „ W (e< - 1)*. - ^L^l) Mv . (27) 

K f Jo Jo v n J 

in (1261) and (1271) . p(a) is given by 

p(a) = - 1) (a',„(^) 2/ " + A^)^/*). (28) 

The proof of Theorem [2] is in Appendix [Bj 

2) Special Case When Interference Experiences General Fading: 

a) General Fading, Noise, a = A: Using the following identity to simplify the results in 
Theorem [2l 




we get the SINR distribution as 



Z f (T) = l- 




(30) 



The mean achievable rate is simplified into 



where 
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b) General Fading, Small Noise, a > 2: As the noise power is approaching zero, we have 
= 1 " P (a)T^R } + ' PM^W ' + ^ (32) 



When a 2 = 0, we further get 



1 _ -p(a)r 2 /«R2 



and the mean achievable rate is 

I f°°l- e -p(«)(e*-i) 2/Q ^ 



3) Special Case When Interference Experiences Rayleigh Fading: For Rayleigh fading, the 
results are essentially the same as that in the general fading case. The only additional simplifi- 
cation is the evaluation of p(a). We just give the result in the very special case when a 2 = 
and a = 4. 

When the interference experiences Rayleigh fading, we have E g (g^) = p Jjj 50 g^e'^dg = 
(1/2) • and consequently, 

p(4) = y(A' m ^^+A'^). (35) 
The SINR distribution is then given by 

Z f (T) = 1 — . (36) 

The mean achievable rate is simplified into 

r°° i- e -p(4)R}Vei=i 

Tf = / o ; dt 

1 Jo p(A)R 2 f y/F^T 

Jo y 2 + p 2 {4)R} 

These expressions are convenient for numerical evaluation. 



C. Mean Achievable Rates of Nonsubscribers and Subscribers 

There are two types of nonsubscribers: those outside nonsubscribers who access MBSs and 
those inside nonsubscribers who access FAPs. When the number of outside nonsubscribers in a 
macrocell is no greater than the total number of subchannels (i.e., U out < M), each nonsubscriber 
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UE exclusively occupies a subchannel, and its mean achievable rate is r m . However, when 
[7 out > M, those U out UEs share the M subchannels with mean achievable rate j¥—T m . Since 

k'out 

the evaluation is conditioned upon the existence of at least one UE, the mean achievable rate of 
an outside nonsubscriber UE is given by 

_ EtU PRut = ^} + YZm + i PRut = ^}f 
Tout l-F{U out = 0} Tm - { j 



Similarly, the mean achievable rate of an inside nonsubscriber is given by 

_ PWn = j} + £°° =Ms+1 nu ia = j}f 



Tf- (39) 



1 - ¥{U- m = 0} 

When averaged over both inside and outside nonsubscribers, the overall mean achievable rate of 
nonsubscriber is obtained as 

XfAinTVJTLj 



^out + AfAinirrCj: 

Regarding subscribers, note that they are exclusively served by FAPs. Similar to the analysis 
of nonsubscribers, the mean achievable rate of a subscriber is given by 

Ts = i-nu s = o} Tf - (41) 

V. Numerical Results 

The numerical results are obtained according to both the analytical results we have derived 
and Monte Carlo simulation. The default configurations of system model are as follows (also see 
Table IJ). The total number of subchannels is M = 20 and the coverage radius of each femtocell 
is Rf = 10m. The transmit power of an FAP is set as Pf = 13dBm and that of an MBS 
P m = 39dBm. We set the path loss exponent as a = 4, with all links experiencing Rayleigh 
fading of normalized /i = 1. The wall penetration loss is set as W = — 6dB. We focus on the 
interference-limited regime, and for simplicity we ignore the noise power (i.e., o 1 = 0). The 
density of MBSs is set as X m = 0.00001 and of FAPs Xf = 0.0001. So the average coverage 
area of an MBS is roughly equal to a circle with a radius of 180m, and on average there are ten 
FAPs within the coverage area of an MBS. Unless otherwise specified, the subscribers and inside 
nonsubscribers are distributed within an FAP coverage are with intensities X s = X m = 0.015. 
The intensity of outside nonsubscribers is set as A out = 0.0001. So on average, each femtocell 
contains about 10 subscribers and inside nonsubscribers, and each macrocell contains about 
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10 outside nonsubscribers. Note that due to statistical fluctuations, the number of UEs may 
occasionally exceeds the number of available subchannels. 

Figure |4] displays the SINR distributions of macrocell UEs and femtocell UEs, when the 
number of shared subchannels is set as M s = 10. We plot in dotted curves the analytical results, 
and we also plot the empirical cdfs obtained from Monte Carlo simulation. The curves reveal 
that the simulation results match the analytical results well, thus corroborating the accuracy of 
our theoretical analysis. From the SINR distributions, we observe that while the macrocell UEs 
experience a fair amount of interference, due to the shrinking cell size, the interference for 
femtocell UEs is substantially alleviated. 

Figure [5] displays the mean achievable rates of macrocell UEs and femtocell UEs. In Figure 
5(a)| the performance is shown as the FAP intensity A/ increases. We observe that the mean 



achievable rates of both macrocell and femtocell UEs only mildly decrease with the increasing 
of A/, suggesting that the network may robustly tolerate the deployment of FAPs. In Figure [5(b)] 
the performance is shown as the outside nonsubscribers intensity A out increases. We observe 
that the mean achievable rates of macrocell and femtocell UEs drop initially and then tend to 
be stable. To interpret this behavior, we note that as the intensity of outside nonsubscribers 
begins to increase, more subchannels become occupied by MBSs, incurring more macrocell tier 
interference; however, when the intensity of outside nonsubscribers is sufficiently large, almost 
all the subchannels are persistently occupied by MBSs with the UEs served by time-sharing, 
and then the interference saturates thus leading to stable performance for UEs. 

Figure [6] displays the mean achievable rates of macrocell UEs and femtocell UEs under 
changing transmission powers. From Figure |6(a)[ we observe that as the transmission power 



of FAPs increases, the rates of femtocell UEs increase correspondingly, and meanwhile the rates 
of macrocell UEs slightly decrease due to the increased cross-tier interference. The opposite 
trends are shown in Figure |6(b)[ where we observe that as the transmission power of MBSs 
increases, the rates of femtocell UEs deteriorate seriously, while the rates of macrocell UEs only 
marginally increase since they are already stuck in the interference-limited regime. 

Figure [7] displays the mean achievable rates of nonsubscribers and subscribers under changing 
number of shared subchannels, in which the rate of nonsubscribers is averaged over both 
outside and inside nonsubscribers. When too few subchannels are to be shared by each FAP 
(i.e., small M s ), the mean achievable rate of nonsubscribers is rather low while subscribers 
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enjoy a good spectral efficiency. At the other extreme, when most of the subchannels are to 
be shared to nonsubscribers, the mean achievable rate of subscribers also deteriorates seriously. 
Nevertheless, from the figure we observe the existence of a stable compromise at which the rates 
of both subscribers and nonsubscribers do not drop much from their maxima. For the default 
configuration in our numerical study, the stable compromise occurs when the value of M s lies in 
the range of [7, 12], corresponding to a reasonably wide tuning range for system designer when 
provisioning the resource. 

Figure [8] displays the mean achievable rate of nonsubscribers and subscribers under changing 
proportions of inside nonsubscribers and subscribers in femtocells. For a fair comparison, we fix 



the sum intensity of the two types of femtocell UEs, as X s + X [n = 0.03. Figure [8(a)] gives the 
performance in the case when most of the femtocell UEs are nonsubscribers, while Figure |8(b)| 
corresponds to the opponent case. From the curves, we observe that in order to achieve a good 
performance for both types of UEs, the number of shared subchannels M s should be adjusted 
based on the intensity of inside nonsubscribers. 

VI. Conclusions 

In this paper, we explored the application of stochastic geometry in the analysis of hybrid 
access mechanisms for multi-channel two-tier networks, focusing on the evaluation of the tradeoff 
between nonsubscribers and subscribers. We characterized several key statistics of UEs and 
subchannels, and established SINR distribution and mean achievable rate for each type of UEs. 

Our analysis revealed the interaction among the various parameters in the network model, and 
thus shed useful insights into the choice of network parameters and the provisioning of resource 
in system design. From our numerical study, we observe that although there is an apparent conflict 
between the interests of nonsubscribers and subscribers, there usually exists a reasonably wide 
tuning range over which nonsubscribers and subscribers attain a stable compromise at which the 
rates of both subscribers and nonsubscribers do not drop much from their maxima. 
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Appendix A 

Let r be the distance between an outside nonsubscriber and its serving MBS. Since an outside 
nonsubscriber always chooses its nearest MBS to access, the cdf of r is obtained by 

F{r <R} = 1 - P{no MBS closer than R} 

= l-e- XmnR2 . (42) 

Then, the pdf of r is given by 

f(r) = e- Xm7rr ' 2 27rX m r. (43) 

Having the pdf of r, we derive the cdf of SINR and the mean achievable rate for an outside 
nonsubscriber UE. Assuming that the considered UE is at distance r from its serving MBS, g. 
is the fading of an interference link, and R. is the distance between the UE and an interfering 
access point, the SINR experienced by the UE is given by SINR = r + J 2 » where I c = 
X^e$' P m giR^ a is the interference from the macrocell tier (excluding the serving MBS 
itself which is denoted by b ) and If = Xlje^ WPfQjRj a is the interference from the femtocell 
tier with the wall penetration loss taken into account. 

We can thus derive the cdf of the SINR following 

Z m {T) = 1-P{SINR>T} 

= 1 - [InX^e-^F { - - > A dr 

Jo {I c + I f + a 2 J 

f°° o ( Tr a ^ 

= 1 - / 27rA m re-" Amr PU>— (/ c + // + ^ 2 ) dr 

JO I °m J 

= 1 -J 27rA m re-^ 2 E(exp (-^(J c + // + a 2 )^j )dr 

= l-/27rA m re-^ 2 -^^ 2 /^/: 7c+J/ (^)dr. (44) 
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By the definition of Laplace transform we get 

4tc+/» = ^ c +/,( e ~ S(/c+//) ) 

= £^(exp(- S ^ P m g r R7^Ev f (eM-sJ2 WP f9j R r< 

= E„ m ( H E g XeM-s9iR7 a P m )))E^]jE gj (eM-mRr WP f)) 
= E* m ( J] L g {sR- a P m ))E^\{L g {sR- a WP f ) 

= exp(-27rA^^ (1 - L g {sv~ a P m ))vdv ) x 



exp ( - 2ttA / / J (l-Lg(sw- a WP f ))wdw^ (45) 

(B) 

where the last equality follows from the pgf of PPP ITT3T1 . Next let us evaluate (A) and (B) 
respectively. 



noo 
(1 - e~ sv - ap ^)f(g)vdgdv 

nOO , POO 



= ~ e ~ SV ' aPm9 > dv )^9)dg. (46) 

Let y = v~ a , then v = y~^, so we get 

i poo pr~ a 

{A) = a I (/ ( 1 - e_aiWl ')v" l " 1 dl/)/U)dp. 
Integrating by parts, since y~*~ 1 dy = —^dy~i, we obtain 

(A) = - ((l-e- sP ^y)y- 2 «{~) r -/ sPmye-^y-^C-^dy)/^)^ 
a J \ 2 o J 2 / 

We have assumed a > 2, so lim^ (l — e~ sPm9y )y~* = 0. Letting t = sP m gy, we obtain 

1 /-oo psP m gr- a 

(A) = -J (-(l-e- sP ^ r ~ a )r 2 + (sP m g^ J e-H l ~- 1 dy)f{g)dg 

r 2 r 2 1 a f°° i / 2 2 \ 

= -- + -^(^ m r- Q ) + -( s p m )" / ^(r(i--)-r(i-- sP^r-"))/^. 

II 2 Jo \ a a / 
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From the properties of Gamma function, we have 



r 











H 




- 2 r( 








a ' 


v a J 



r ( 1 - sP^r"^ = ( S P^r- Q )- 2 / a e- sP '^" a - -r ( --, sP m gr- 
\ a J a \ a 

so that (A) can be simplified as 

(A) = -^ + ^P m )^M(n--,sP m gr- a )-T(--))). 
2 a V V a a / / 

The evaluation of (5) is almost the same, and we have 

(5) = --(sWP f ^T(--)E 9 (g^. 

a a 

Substituting (A) and (B) into £/ o+ /„(s) and rearranging terms in Z m (T) hence leads to ©. 

Now we evaluate the mean achievable rate. Since for a positive random variable X, E[X] = 
J P{X > t}dt, we have 



P[ln(l + SINR)] 



,-?rA m r 2 



In 1 + 



OO /"00 

„-7rA m r 2 



I c + I f + ° 2 
Pmhr~ 



27rAr„rdr 



> t >dt2ir\ m rdr 



OO /"OO 

„-7rA m r 2 



OO /"OO 



OO /"OO 



JO 



™ f , / r m Air 
P{ ln(l + — jp- 
L i c + if 

/ F\h > —(I c +I f + (7 2 )(e* - 1) [dt27rA m rdr 

70 Pyn ' 

E(e^^ (Jc+ ^ +ff2)(e '- 1) )dt2 7 rA m rdr 
^(e*-i), /^r a (e*-l)- 



dt27rA rr ,rdr. 



(47) 



From the derivation of Z m {T) we have already obtained the expression of £>i c +i f 
When substituting into (l47l) . we arrive at (fTOl) and thus establish Theorem [IJ 



Appendix B 



Let r be the distance between a femtocell UE and its serving FAR Because the femtocell UEs 
are uniformly distributed in the circular coverage area of radius Rf of each FAP, the pdf of r is 
given by /(r) = 2r/R'j. 
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Denote by I' c and I'j the interference strengths from the macrocell tier and the femtocell tier 
respectively. Similar to the proof of Theorem [Q we have 

W e~ Pf £ ^(^)dr. (48) 

Following the same procedure as in the derivation of (1431 ), we get the Laplace transform of 

I' c + I'j as 

, />CO POO 

= exp ( - 2vrA^ y (1 - £ fl (siT'W.P m ))T;di; - 2ttAJ y (1 - L g (sw- a W 2 P f ))wdw 

= exp ( - 7rsSr(l - ^)(AU^P m )« + A'^P/^E^)). (49) 

In the above, it is noteworthy that the interference from an FAP penetrates two walls thus the 
loss becoming W 2 instead of W. Substituting ( |49b into (|48l) with v = r 2 , we get the SINR 
distribution in Theorem El Similar to (|47l) . the mean achievable rate is given by 

= ^ ^ e / L Pc+I , ( p/ J ^dtdr. (50) 

Substituting (|49l ) into (l50l) . we get the mean achievable rate in Theorem |2l 
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TABLE I 
SYSTEM PARAMETERS 



Symbol 


Description 


Value 




PPPs defining the MBSs and FAPs 


N/A 




Density of MBSs 


0.00001 MBS/m 2 


A/ 


Density of FAPs 


0.0001 FAP/m 2 


«/ 


Radius of femtocell 


10m 


Aout 


Densities of outside nonsubscribers 


0.0001 user/m 2 


A s , Ai n 


Densities of subscribers and inside nonsubscribers 


0.015 user/m 2 




Transmit power at MBS and FAP 


39dBm,13dBm 


M 


Number of subchannels at each access point 


20 


M s , M r 


Number of subchannels shared and reserved by each femtocell 


Not fixed 


a 


Path loss exponent 


4 


P- 


Rayleigh fading parameter 


1 (normalized) 


W 


Wall penetration loss 


-6dB 


a 2 


Noise power 


(interference-limited regime) 


pm pf 
-'busy 5 -'busy 


Probabilities that a given subchannel is used by an MBS and an FAP 


Not fixed 




PPPs defining the MBSs and FAPs that interfere a given subchannel 


N/A 


Ami A'y 


Densities of MBSs and FAPs that interfere a given subchannel 


Not fixed 


Uout j CAn 


Numbers of nonsubscribers that access a given MBS and a given FAP 


Not fixed 


17. 


Numbers of subscribers that access a given FAP 


Not fixed 


T/, T m 


Mean achievable rates of femtocell UEs and macrocell UEs 


Not fixed 




Mean achievable rates of subscribers and nonsubscribers 


Not fixed 
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Fig. 1. Downlink two-tier network model for hybrid access femtocell. 



Network Model 




Fig. 2. The Voronoi macrocell topology, in which each Voronoi cell is the coverage area of a macrocell and each small circle 
represents a femtocell coverage area. 
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subchannels 



FAP 1 



FAP 2 



FAP 3 



Subchannels allocated to the nonsubscribers 
Subchannels reserved for the subscribers 



Fig. 3. Spectrum allocation in each hybrid access femtocell. All the M s shared subchannels are randomly selected by each 
FAP. 



SINR Distribution 




SINR (dB) 



Fig. 4. Cdfs of SINRs for macrocell UEs and femtocell UEs, when M s — 10. 
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16 - 
14 



- Simulative Macrocell UEs 

- Theoretical Macrocell UEs 

- Simulative Femtocell UEs 
Theoretical Femtocell UEs 



2 3 
Density of the FAPs X f 



- Simulative Macrocell UEs 

- Theoretical Macrocell UEs 

- Simulative Femtocell UEs 
Theoretical Femtocell UEs 



S 1" - 




1 2 3 

The density ot the outside nonsubscribers X 



(a) Change the intensity of FAPs A/. (b) Change the intensity of outside nonsubscribers A ou t- 

Fig. 5. Mean achievable rates of macrocell UEs and femtocell UEs with changing intensities of FAPs and outside nonsubscribers. 




The transmission power of the FAPs (dBm) The transmission power of the MBSs (dbm) 



(a) Change Pf, P m = 39dBm fixed. (b) Change P m , P f = 13dBm fixed. 

Fig. 6. Mean achievable rates of macrocell UEs and femtocell UEs with changing transmission powers. 
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Fig. 7. Performance of nonsubscribers and subscribers with changing number of shared subchannals M s . 
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Multi-channel Hybrid Access Femtocells: 
A Stochastic Geometric Analysis 

Yi Zhong and Wenyi Zhang, Senior Member, IEEE 



Abstract 

[TO UPDATE] For femtocells, the channel access mechanism can be configured to be open access 
or closed access. Open access may lead to performance deterioration for femtocell subscribers, while 
closed access may cause serious interference to nonsubscribers served by macrocells. A compromise 
solution is hybrid access, in which only a fraction of the available resource is shared to nonsubscribers 
while the remaining reserved for subscribers. This paper focuses on a hybrid access mechanism for multi- 
channel femtocells which employ orthogonal spectrum access schemes, like OFDMA. Considering a 
randomized channel assignment strategy, we analyze the performance in the downlink. Using stochastic 
geometry as technical tools, we derive the blocking probabilities, mean achievable rates, and distributions 
of the signal-to-interference-plus-noise ratio of both nonsubscribers and subscribers. For these key 
performance indicators, we establish either explicit or integral closed form expressions, which are 
amenable to numerical evaluations. The analysis reveals the performance tradeoff between subscribers 
and nonsubscribers, and is corroborated by numerical simulations. 

Index Terms 

Channel management, femtocell, hybrid access, spatial Poisson process, two-scale approximation, 
two-tier network 
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1 

I. Introduction 

In current cellular network services, about 50% of phone calls and 70% of data services 
take place indoors [?]. For such indoor use cases, network coverage is a critical issue. One 
way to improve the indoor performance is to deploy the so-called femtocell access points 
(FAPs), also known as home base stations, besides macrocell base stations (MBSs). Femtocells 
are small cellular base stations, typically designed for use in home or small business [?][?]. 
The use of femtocells not only benefits the users, but also the operators. As the distance 
between the transmitter and receiver is reduced, users will enjoy high quality links and power 
savings. Furthermore, the reduced transmission range also creates more spatial reuse and reduces 
electromagnetic interference. 

Among the many challenges faced by femtocells, and more generally, two-tier networks, is 
the issue of interference; see Figure [Q The two-tier interference problem differs from that in 
traditional single-tier networks in several important aspects: First, due to the limitations of access 
mechanism, a user equipment (UE) may not be able to connect to the access point which offers 
the best service. Second, since femtocells connect to operator's core network via subscribers' 
private ISP, coordination between macrocells and femtocells and among femtocells is limited. 
Finally, compared to planned macrocell deployments, femtocells are usually deployed in an ad 
hoc manner, and the thus randomly placed femtocells make it difficult to manage the interference. 
In two-tier networks comprising of macrocells and femtocells, interference can be categorized 
into two types: 

• Cross-tier: This refers to the interference from one tier to the other; that is, femtocell 
interference at a macrocell and macrocell interference at a femtocell. 

• Co-tier: This refers to the interference within a tier; that is, femtocell interference at a 
nearby femtocell and macrocell interference at a nearby macrocell. 

In this paper, we consider two-tier networks that are based on multicarrier techniques, for 
example those deploying LTE or WiMAX standards, which use orthogonal frequency-division 
multiple access (OFDMA) techniques. In multicarrier systems, the available spectrum is divided 
into orthogonal subcarriers, which are then grouped into multiple subchannels, assigned to dif- 
ferent users. Due to the flexibility in channel assignment, both cross-tier and co-tier interferences 
may be alleviated. 
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The access mechanism of femtocells (see, e.g., [?]) is a key factor that affects the perfor- 
mance of two-tier networks, and generally can be classified as follows, where we call the UEs 
that register to a femtocell as subscribers, and those that do not register to any femtocell as 
nonsubscribers. 

• Closed access: An FAP only allows its subscribers to connect. 

• Open access: An FAP allows all its covered UEs, no matter registered or not, to connect. 

• Hybrid access: An FAP allows its covered nonsubscribers to connect via a subset of its 
available subchannels, and reserves the remaining subchannels for its subscribers. 

Hybrid access [?] is an intermediate approach, in which a fraction of the total available resource 
is allocated to nonsubscribers. By doing so, nonsubscribers near an FAP may handover into the 
femtocell to avoid high interference; meanwhile, with a certain amount of resource reserved 
for subscribers, the performance of subscribers may be well assured even in the presence of 
nonsubscribers. 

In hybrid access, a central issue is how to allocate the resource between subscribers and 
nonsubscribers. Previous studies [?] [?] indicate that hybrid access improves the network perfor- 
mance at the cost of reduced performance for subscribers, therefore suggesting a tradeoff between 
the performance of nonsubscribers and subscribers. In this paper, we consider a hybrid access 
mechanism that uses a randomized channel assignment strategy, and analyze the performance in 
the downlink of both macrocells and femtocells. We employ stochastic geometry to characterize 
the spatial distributions of users as well as access points; see, e.g., [?] and references therein for 
its recent applications in wireless networks. Accordingly, we derive the key performance indi- 
cators including blocking probabilities, mean achievable rates, and distributions of the signal-to- 
interference-plus-noise ratios (SINRs), of both nonsubscribers and subscribers. In our study, we 
establish general integral expressions for the performance indicators, and closed form expressions 
under specific model parameters. With the obtained results, we reveal how the performance of 
subscribers and nonsubscribers trades off each other. 

An overview of related works is as follows. In [?], the authors proposed to study key per- 
formance indicators for cellular networks, such as coverage probabilities and mean achievable 
rates, through models and tools in stochastic geometry. Our work explores that approach in hybrid 
access femtocell networks, focusing on multi-channel models and the downlink performance. In 
[?], the considered scheme divides the spectrum resource into two orthogonal parts which are 
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assigned to macrocells and femtocells, respectively, with femtocells being closed access. In [?], 
the authors considered two-tier femtocell networks using time-hopped CDMA, examining the 
uplink outage probability and the interference avoidance capability. In [?] and [?], the authors 
studied the performance of various femtocell access mechanisms, under substantially different 
system models from ours. Our work examines the interaction among multiple MBSs and FAPs, 
and focuses on scenarios in which macrocells and femtocells share the same frequency band, 
following a hybrid access mechanism, as will be elaborated in Section UH 

The remaining part of the paper is organized as follows. Section [n] describes the two-tier 
network model, the channel assignment strategy, and the hybrid access mechanism. Based on 
a two-scale approximation for the spatial distributions of femtocells, Section [III] analyzes the 
statistical behavior of UEs, deriving the distributions of the number of UEs connecting to either 
an FAP or an MBS, as well as the probabilities of a subchannel being used by either an FAP 
or an MBS. Built upon those statistics, Section [IV] establishes expressions for the distributions 
of the SINRs, and the mean achievable rates. Section |V] illustrates the aforementioned analysis 
by numerical results, which are also corroborated by simulations. Finally, Section IVIl concludes 
the paper. 

II. Network Model 

A. Hybrid Access Femtocells 

In the two-tier network, we consider two types of access points, MBSs and FAPs. The MBSs 
constitute the macrocell tier, and they induce a Voronoi tessellation of the plane (see Figure [2]). 
When a UE attempts to access the macrocell network, it chooses to connect to the MBS in the 
Voronoi cell which it is situated in. An FAP provides network access to UEs in its vicinity, and 
we assume that all FAPs have a covering radius of Rf. Within the covered circular area of each 
FAP are two types of UEs, called subscribers and inside nonsubscribers. Inside nonsubscribers 
are those UEs who gather around an FAP without subscribing to its service; for example, transient 
customers in a shop or restaurant. Besides those two types of UEs, we also consider a third type 
of UEs, outside nonsubscribers, who are uniformly scattered over the whole plane, corresponding 
to those regular macrocell network users. 

The available spectrum is evenly divided into M subchannels, which are to be shared by both 
the macrocell tier and the femtocell tier. Each FAP is configured to allocate a fixed number, M s , 
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of subchannels for its covered inside nonsubscribers. These M s subchannels are called shared 
subchannels, and the remaining M r = M — M s subchannels are called reserved subchannels 
as they are reserved for the subscribers. In the considered hybrid access mechanism, each FAP 
selects its shared subchannels randomly, and independently of other FAPs. We assume that 
each UE, whether subscriber or nonsubscriber, needs one subchannel for accessing the two-tier 
network. When a UE accessing an MBS or an FAP, the serving subchannel is selected randomly 
(see Figure [3]). 

The hybrid access mechanism operates as follows. 

• A subscriber accesses to one of the M r reserved subchannels of its corresponding FAP. 
When there are more than M r subscribers in an FAP, they are served by time-sharing with 
equal time proportion. 

• An inside nonsubscriber accesses to one of the M s shared subchannels of its covering 
FAP. When there are more than M s inside nonsubscribers in an FAP, they are served by 
time-sharing with equal time proportion. 

• An outside nonsubscriber accesses the MBS in the Voronoi cell which the outside nonsub- 
scriber is situated in. When there are more than M outside nonsubscribers in the Voronoi 
cell, they are served by time-sharing with equal proportion. 

B. Mathematical Model and Two-scale Approximation 

To formulate the aforementioned hybrid access scenario mathematically, we model the spatial 
distributions of the nodes using spatial Poisson processes as follows. The MBSs constitute a 
homogeneous Poisson point process (PPP) $ m of intensity A m on the plane; the FAPs constitute 
another homogeneous PPP $j of intensity A/. In the circular covered area of radius Rt of each 
FAP, the subscribers are distributed according to a homogeneous PPP of intensity A s , and the 
inside nonsubscribers are distributed according to another homogeneous PPP of intensity A in . 
Outside nonsubscribers constitute on the whole plane a homogeneous PPP of intensity A out . All 
the PPPs are mutually independent. 

In this paper, we focus on the downlink performance. The transmit power is set to a constant 
value P m for an MBS, and Pf for an FAP. For the sake of convenience, we adopt a standard path 
loss propagation model with path loss exponent a > 2. Regarding fading, we assume that the 
link between the serving access point (either an MBS or an FAP) and the served UE experiences 
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Rayleigh fading with parameter fi. The received signal power of a UE at a distance r from its 
serving access point therefore is P m hr~ a (MBS) or Pfhr~ a (FAP) where h ~ Exp(/i). The 
fading of interference links may follow an arbitrary probability distribution, and is denoted by 
g. Furthermore, considering the typical scenario of indoor femtocell deployment, we introduce a 
wall isolation at the boundary of each FAP coverage area, which incurs a wall penetration loss 
factor W < 1. For all receivers, the noise power is a 2 . 

The different PPPs corresponding to different entities in the network interact in a complicated 
way, thus making a rigorous statistical analysis extremely difficult. For example, an inside 
nonsubscriber may be covered by more than one FAPs, thus leading to the delicate issue of 
FAP selection, and furthermore rendering the subchannel usage distributions among FAPs and 
MBSs intrinsically correlated. To overcome the technical difficulties due to spatial interactions, in 
the subsequent analysis we propose a two-scale approximation for the network model, motivated 
by the fact that the covered area of an FAP is significantly smaller than that of an MBS. The 
two-scale approximation consists of two views, the macro-scale view and the micro-scale view. 
The macro-scale view concerns an observer outside the coverage area of an FAP, and in that 
view the whole coverage area of the FAP shrinks to a single point, marked by the numbers 
of subscribers and inside nonsubscribers therein. The micro-scale view concerns an observer 
inside the coverage area of an FAP, and in that view the coverage area is still circular with 
radius Rf in which the subscribers and inside nonsubscribers are spatially distributed. By such 
a two-scale approximation, an inside nonsubscriber can only be covered by a unique FAP, and 
the coverage area of an FAP can only be within a unique Voronoi cell of an MBS. These 
consequences substantially simplify the performance analysis. In Section |V] we validate the 
two-scale approximation method through comparing analytical results and simulation results, for 
network parameters of practical interest. 

III. Statistics of UEs and Subchannels 

In this section, we characterize the distributions of UEs connecting to different types of access 
points, and the distributions of used subchannels in MBSs and FAPs. The analysis is based on a 
snapshot of the network model, and the obtained results will then be applied for characterizing 
the distributions of SINRs and achievable rates in Section HVl 
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A. Distributions of UEs 

Let U s be the number of subscribers accessing a given FAP and from our model we have 
U s ~ Poisson(A s 7ri?j). Similarly, let U m be the number of inside nonsubscribers accessing a 
given FAP, and we have U in ~ Poisson(A in 7ri?j). 

The number of outside nonsubscribers who access a given MBS, denoted by U out , is character- 
ized as follows. We note that the macrocell coverage area is a Voronoi cell, and denote by S the 
area of the Voronoi cell. There is no known closed form expression of the probability density 
function (pdf) of S, whereas a simple approximation [?] has proven sufficiently accurate for 
practical purposes. Considering scaling, the approximate pdf of the size of a macrocell coverage 
area is given by 

f(S) = ^/X(SA m )§ exp(-^A m )A m . (1) 

Conditioned upon S, the number of outside nonsubscribers is a Poisson random variable with 
mean A^tS*. The probability generating function (pgf) of the unconditioned U out is thus given 

by 

G(z) = ^°°exp (A out (z - l)s)f(S)dS. (2) 

Plugging in the approximate pdf of S and simplifying the integral, we get 

343 /7/7 A out yl 

The distribution of U out is therefore given by the derivatives of G(z), 

, GW(0) 

P{^out=0 = — r 1 , i = 0,l,.... (4) 

B. Distributions of Subchannel Usage 

Since the subchannels are uniformly and independently selected by each FAP, it suffices to 
analyze an arbitrary one of them. Let us examine the probability that a given subchannel is used 
by an MBS or an FAP. In order to obtain the probability, we evaluate the average number of 
subchannels used by an MBS or an FAP, and then normalize the average number by the total 
number of subchannels, M. 
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The probability that a subchannel is used by an FAP is 

^ oo oo 

Pbusy,/ = ^(Yl mil1 ^' M rM U * =0 + X] mil1 ^' M sMUin = j} 
i=0 j=0 

.M r M OO 

= = o + a - if) E p {^ = o + 



«=0 i=M r +l 
-.M s oo 

j=0 j=M s + l 

For a Poisson random variable iV ~ Poisson(A), its cumulative distribution function (cdf) is 
given by 

± m = (i = ±^ = *±hx>, 

i=0 i=0 

where T(s,x) = f°° i s_1 e~*dt is the incomplete gamma function. Using © to simplify -Ptmsy,/* 
we get 



Pbus y ,/ = 1 + (A s 7r J Rjr(M r , A^Pj) - T(M r + 1, A.ttPJ)) + 

+ ^T4T(A i „7r J Rjr(M s ,A in 7ri?J) -r(M s + l,A in 7ri?J)). (7) 



M M s 

The probability that a subchannel is used by an MBS is 

j oo 

n usy , m = — E min {^^}p{^out = ^}, (8) 

i=0 

where P{C/ out = z} is given by ©. 

The spatial process of FAPs that use a given subchannel is the independent thinning of the 
original PPP of FAPs $j by the probability -Pbusy,/* denoted by with intensity A^ = XfPbusyj- 
The term "independent thinning" means that $^ can be viewed as obtained from by inde- 
pendently removing points with probability 1 — Pbusy,/- Similarly, the spatial process of MBSs 
that use a given subchannel is the independent thinning of the original PPP of MBSs $ m by the 
probability .Pbusy,™, denoted by <3>' m with intensity X' m = A m Pbu Sy ,m- These two independently 
thinned PPPs will prove useful in the subsequent analysis. 

IV. Analysis of SINR and Mean Achievable Rate 

In this section, we derive the distributions of SINRs and the mean achievable rates of all the 
three types of UEs. For each type of UEs, we begin with general settings, and then simplify the 
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general results under several specific parameters to gain insights. The mean achievable rates are 
the averaged instantaneous achievable rates over both channel fading and spatial distributions of 
UEs and access points. 

A. Macrocell UEs 

1) General Case: For an active UE served by an MBS, it must be occupying one subchannel 
of the MBS. The following theorem gives the cdf of SINR and the mean achievable rate of each 
active macrocell UE in general case. 

Theorem 1: The cdf of the SINR of a macrocell UE, denoted by Z m (T) = P{SINR < T}, 
is given by 

Z m (T) = 1 - nxj exp (nv\' m (l - /3(T, a) - -^—) - (9) 

Jo ' busy,m m ' 

The mean achievable rate of a macrocell UE is given by 

exp (nv\' m (l - /3(e* - l,a) - ) - ^ ± } -)dtdv. (10) 

^ busy,m * ra ' 

In © and fl©, (3{T, a) is given by 

ma) = ^ E Li(r ( - 2 -,,T g) - (i + 

The proof of Theorem [T] is in Appendix [A] 

In Theorem [H Z m (T) in © gives the probability that the SINR is below a given target level 
T, and r m in (TTOl) gives the mean achievable rate of a macrocell UE. The integrals in © and 
(flOl) can be evaluated by numerical methods, and furthermore, they can be simplified to concise 
forms in a number of special cases. The cases of interest in this paper are the combinations of 
(i) interference experiencing Rayleigh fading, (ii) the path loss exponent being a = 4, and (iii) 
no noise or the noise is weak, a 2 — > 0. The following discussions are categorized into two parts 
according to the interference fading. 

2) Special Case When Interference Experiences General Fading: 

a) General Fading, Noise, a = 4: We can simplify © via invoking the identity 
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where Q(x) = h= J r ^° exp(— y 2 /2)dy is the standard Gaussian tail probability. Letting a(T) 
7rA' m (/3(T,4) + X m /X' m - 1), we have 



ZJV = 1 - ,A ml /^exp (^) Q ( a(r) ^). (13) 

For a given distribution of interference fading, j3(T, 4) can be evaluated by numeric methods, 
and subsequently the SINR distribution Z m (T) can be obtained conveniently. 

b) General Fading, Small Noise, a > 2: Consider the interference-limited regime in which 
the noise power is approaching zero. Using the expansion e~ x = 1 — x + o(x) as x — > 0, the 
SINR distribution behaves like 



2\ 



m( } = " 1 + (B(T a) - 1) P h + 7 / xx«/2+i + °( a 

(14) 

When a 2 = 0, the SINR distribution is further simplified into 

Zm ^ = 1 - ,+(r(tJ up • (15) 

1 + {p[l,a) - l)-Pb US y, m 

3) Special Case When Interference Experiences Rayleigh Fading: Here we consider the case 
where the interference experiences Rayleigh fading with mean i.e. g ~ Exp(/i). In this case, 
the results are simplified as follows. 

Corollary 1: When the interference follows Rayleigh fading, the cdf of the SINR is 

Z m (T) = 1 - n\ m J exp ( - irv\ m - irv\' m ip(T, a, T~ 2/a ) - 

fPfWT \ ^V\. 
_7 ™ A M p — p ) dv - ( 16 ) 

The mean achievable rate is 

= ^X m J J exp ( - 7ruA m - irvX^ip^e* - I, a, (e* - l) _2/o ) - 

—7Tv\f(p ( ^ i,a,0j Mdudt, (17) 



7~n 



where 



The results in Corollary Q] are further simplified in the following special cases. 



(18) 
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a) Rayleigh Fading, Noise, a = 4: Using (fT2l) . we get 



Zm{T) = 1 - .A ra ,/^exp (^) Q ( K (T)^). (19) 
The mean achievable rate of a macrocell UE is simplified into 



'/^^^ ex KC^^)^(^'-Vw^y )d ' <20) 

where 

k(T) = Ti\ m + 7rA' m VT(7r/2 - arctan VT) + p ■ (21) 

^ V -*ra 

Z?) Rayleigh Fading, No Noise, a > 2: When a 2 = 0, the SINR distribution is simplified 

into 

Z m (T) = 1 -, . (22) 



The mean achievable rate of a macrocell UE is simplified into 

1 



1 + ^(e* - 1, a, (e* - !)-*/«) + ^>(™?^, a, 0) 



dt. (23) 



Specifically, when a = 4 we obtain 



z ™ ( T ) = 1 7 v TwWS • (24) 

1 + y/T{* - arctan VT + z^*£l)P busy>m 



The mean achievable rate of a macrocell UE is simplified into 

2 



tan?/ + (tt/2 - y)P b u Sy ,m + f jjr- 



:dy. (25) 



From a practical perspective, it is desirable to shape Z m (T) to make it small for small 
values of T. From ([24]) . we see that there are two approaches to shape Z m (T). First, Z m (T) 
decreases as Pbus y ,m, the probability that a subchannel is used by an MBS, decreases. This may 
be interpreted as an effect of frequency reuse. Second, Z m (T) decreases as the whole term, 
y-yW^^ decreases, which corresponds to a number of network parameters, representing the 
effect due to the deployment of the femtocell tier. 
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B. Femtocell UEs 

1 ) General Case: A UE served by an FAP occupies one subchannel of the FAR The following 
theorem gives the cdf of SINR and the mean achievable rate of each active femtocell UE in 
general case. 

Theorem 2: The cdf of the SINR of a femtocell UE in general case is given by 

Z f( T ) = 1 -Jp f l ex P ( - P{oc)T 2/a v - * ^ )dv. (26) 
The mean achievable rate of a femtocell UE is given by 

V - 4 f ? f - P ( - „ W (e< - 1)*. - ^L^l) Mv . (27) 

K f Jo Jo v n J 

in (1261) and (1271) . p(a) is given by 

p(a) = - 1) (a',„(^) 2/ " + A^)^/*). (28) 

The proof of Theorem [2] is in Appendix [Bj 

2) Special Case When Interference Experiences General Fading: 

a) General Fading, Noise, a = A: Using the following identity to simplify the results in 
Theorem [2l 




we get the SINR distribution as 



Z f (T) = l- 




(30) 



The mean achievable rate is simplified into 



where 
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b) General Fading, Small Noise, a > 2: As the noise power is approaching zero, we have 
= 1 " P (a)T^R } + ' PM^W ' + ^ (32) 



When a 2 = 0, we further get 



1 _ -p(a)r 2 /«R2 



and the mean achievable rate is 

I f°°l- e -p(«)(e*-i) 2/Q ^ 



3) Special Case When Interference Experiences Rayleigh Fading: For Rayleigh fading, the 
results are essentially the same as that in the general fading case. The only additional simplifi- 
cation is the evaluation of p(a). We just give the result in the very special case when a 2 = 
and a = 4. 

When the interference experiences Rayleigh fading, we have E g (g^) = p Jjj 50 g^e'^dg = 
(1/2) • and consequently, 

p(4) = y(A' m ^^+A'^). (35) 
The SINR distribution is then given by 

Z f (T) = 1 — . (36) 

The mean achievable rate is simplified into 

r°° i- e -p(4)R}Vei=i 

Tf = / o ; dt 

1 Jo p(A)R 2 f y/F^T 

Jo y 2 + p 2 {4)R} 

These expressions are convenient for numerical evaluation. 



C. Mean Achievable Rates of Nonsubscribers and Subscribers 

There are two types of nonsubscribers: those outside nonsubscribers who access MBSs and 
those inside nonsubscribers who access FAPs. When the number of outside nonsubscribers in a 
macrocell is no greater than the total number of subchannels (i.e., U out < M), each nonsubscriber 
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UE exclusively occupies a subchannel, and its mean achievable rate is r m . However, when 
U out > M, those U out UEs share the M subchannels with mean achievable rate j¥-r m . Since 
the evaluation is conditioned upon the existence of at least one UE, the mean achievable rate of 
an outside nonsubscriber UE is given by 

_ ggi nuou t = ^} + eZm+i nuout = oj 

r ° ut i-p{f/ out = o} Tm - 1 j 

Similarly, the mean achievable rate of an inside nonsubscriber is given by 

Tin = 1 - w{u- m = o} r/ - (39) 

When averaged over both inside and outside nonsubscribers, the overall mean achievable rate of 
nonsubscriber is obtained as 

r ™ = — x — — • ( 4 °) 

A out ~r AfAinlTrtj! 

Regarding subscribers, note that they are exclusively served by FAPs. Similar to the analysis 
of nonsubscribers, the mean achievable rate of a subscriber is given by 

Ts = 1-F{U. = 0} T/ - (41) 

V. Numerical Results 

The numerical results are obtained according to both the analytical results we have derived 
and Monte Carlo simulation. The default configurations of system model are as follows. The 
total number of subchannels is M = 20 and the coverage radius of each femtocell is Rf — 10m. 
The transmit power of an FAP is set as Pf = 13dBm and that of an MBS P m = 39dBm. We 
set the path loss exponent as a = 4, with interference links experiencing Rayleigh fading. The 
wall penetration loss is set as W = — 6dB. We focus on the interference-limited regime, and for 
simplicity we ignore the noise power (i.e., a 2 = 0). The density of MBSs is set as A m = 0.00001 
and of FAPs A/ = 0.0001. So the average coverage area of an MBS is roughly equal to a circle 
with a radius of 180m, and on average there are ten FAPs within the coverage area of an MBS. 
Unless otherwise specified, the subscribers and inside nonsubscribers are distributed within an 
FAP coverage are with intensities \ s = A in = 0.015. The intensity of outside nonsubscribers is 
set as A out = 0.0001. So on average, each femtocell contains about 10 subscribers and inside 
nonsubscribers, and each macrocell contains about 10 outside nonsubscribers. Note that due to 
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statistical fluctuations, the number of UEs may occasionally exceeds the number of available 
subchannels. 

Figure |4] displays the SINR distributions of macrocell UEs and femtocell UEs, when the 
number of shared subchannels is set as M s = 10. We plot in dotted curves the analytical results, 
and we also plot the empirical cdfs obtained from Monte Carlo simulation. The curves reveal 
that the simulation results match the analytical results well, thus corroborating the accuracy of 
our theoretical analysis. From the SINR distributions, we observe that while the macrocell UEs 
experience a fair amount of interference, due to the shrinking cell size, the interference for 
femtocell UEs is substantially alleviated. 

Figure [5] displays the mean achievable rates of macrocell UEs and femtocell UEs. In Figure 
5(a)| the performance is shown as the FAP intensity A/ increases. We observe that the mean 



achievable rates of both macrocell and femtocell UEs only mildly decrease with the increasing 
of A/, suggesting that the network may robustly tolerate the deployment of FAPs. In Figure [5(b)] 
the performance is shown as the outside nonsubscribers intensity A out increases. We observe 
that the mean achievable rates of macrocell and femtocell UEs drop initially and then tend 
to be stable. This is because the rise of the number of outside nonsubscribers results in more 
subchannels in use at each MBS, which brings more interference to the network. However, when 
the density of outside nonsubscribers is large enough, almost all the subchannels of the MBSs 
are used and the UEs are served by time-sharing with equal time proportion. In this case, the 
interference will not increase any more with the rise of the density of outside nonsubscribers. 
Therefore, for an active UE, the mean rate tends to be stable. 

Figure [6] gives the mean achievable rates of macrocell UEs and femtocell UEs under various 



transmission powers. From Figure |6(a)[ we see that as the transmission power of the FAPs 
increases, the rate of femtocell UEs increases correspondingly. However, the increasing power 
of the FAPs will impose more interference to the macrocell UEs and degrade the performance 



of them. Figure |6(b)| shows that when fixing the power of FAPs and increasing the power of 



MBSs, the performance of femtocell UEs deteriorates seriously. 

Figure |7] gives the performance of nonsubscribers and subscribers under various numbers 
of shared subchannels in an FAP. The mean achievable rate of nonsubscribers is the average 
over both outside and inside nonsubscribers. When M s is small, which is the case that few 
subchannels are shared by each FAP, the mean rate of nonsubscribers is terrible although the 



August 8, 2011 



DRAFT 



15 



performance of subscribers is pretty good. As the number of shared subchannels increases, the 
performance of nonsubscribers improves and of subscribers decreases, which is consistent with 
the actual. At the extreme case that almost all subchannels are shared to the nonsubscribers, the 
performance of subscribers deteriorates severely. This picture gives an evidence of the conflict 
between performance of subscribers and nonsubscribers. A compromise must be taken between 
the performance of both kinds of UEs. According to Figure [/J at the default configuration, it is 
advisable to choose the value of M s in the range of [7, 12]. The result indicates that within a 
large range of the value M s , the performance is pretty good. 

Figure [8] shows the mean achievable rate of nonsubscribers and subscribers under different 
proportions of inside nonsubscribers and subscribers in femtocells. The sum density of all kinds 



of UEs in a femtocell is fixed, which is A s + A in = 0.03. Figure [8(a)] gives the performance in the 
case when most of the femtocell UEs are nonsubscribers while Figure |8(b)| gives the opponent 
case. The results indicate that in order to achieve a good performance for both kinds of UEs, 
the shared number of subchannels in each FAP should changed with the density of the inside 
nonsubscribers. 



VI. Conclusions 

This paper has proposed a tractable network model to evaluate the tradeoff between the 
nonsubscribers and subscribers. We derived the closed form expressions of the SINR distribution 
and the mean achievable rate for each UE. In the special cases, we saw that the expressions 
can be simplified into some quite concise forms. With these closed form results, the relationship 
between the performance of each US and the proportion of shared resource became clear. 

The numerical results show that there is an apparent conflict between the performance of 
nonsubscribers and that of subscribers. The increase of the shared resource will lead to the im- 
prove of the performance of nonsubscribers, but the performance of subscribers will deteriorate. 
Moreover, it is shown that the amount of shared spectrum has a large value space in order to 
balance the performance of both nonsubscribers and subscribers. 
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Appendix A 

Let r be the distance between an outside nonsubscriber and its serving MBS. Since an outside 
nonsubscriber always chooses its nearest MBS to access, the cdf of r is obtained by 

F{r <R} = 1 - P{no MBS closer than R} 

= l-e- XmnR2 . (42) 

Then, the pdf of r is given by 

f(r) = e- Xm7rr ' 2 27rX m r. (43) 

Having the pdf of r, we derive the cdf of SINR and the mean achievable rate for an outside 
nonsubscriber UE. Assuming that the considered UE is at distance r from its serving MBS, g. 
is the fading of an interference link, and R. is the distance between the UE and an interfering 
access point, the SINR experienced by the UE is given by SINR = r + J 2 » where I c = 
X^e$' P m giR^ a is the interference from the macrocell tier (excluding the serving MBS 
itself which is denoted by b ) and If = Xlje^ WPfQjRj a is the interference from the femtocell 
tier with the wall penetration loss taken into account. 

We can thus derive the cdf of the SINR following 

Z m {T) = 1-P{SINR>T} 

= 1 - [InX^e-^F { - - > A dr 

Jo {I c + I f + a 2 J 

f°° o ( Tr a ^ 

= 1 - / 27rA m re-" Amr PU>— (/ c + // + ^ 2 ) dr 

JO I °m J 

= 1 -J 27rA m re-^ 2 E(exp (-^(J c + // + a 2 )^j )dr 

= l-/27rA m re-^ 2 -^^ 2 /^/: 7c+J/ (^)dr. (44) 
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By the definition of Laplace transform we get 

= ^(exp(-s ]T P m ^ J R- Q ))^ / (exp(- S J] WP f9j Rj a j) 

= ^( n ^xex P (-^^rnn)))^(n^( ex p(-^ i? 7 a ^ p /)) 
= E„ m ( n £ g ( SJ Rrp m ))^(n^(^7 Qwp /)) 

= exp(-27rA^^ (1 - £ 9 (stT a P m ))t;dt; ) x 



exp ( - 2vrA / / J (l-L g (sw- a WP f ))wdw^), (45) 

(B) 

where the last equality follows from the pgf of PPP [?]. Next let us evaluate (A) and (B) 
respectively. 

poo poo 

(A) = / (l-e-™- ap «*)f(g)vdgdv 

Jr JO 

= J" {J°°0--e- sv ~ aPm9 )vdv > jf(g)dg. (46) 
Let y = v~ a , then v = so we get 

1 poo f v a 

(A) = -J (J (l-e- sP ^y)y-^- 1 dy)f(g)dg. 
Integrating by parts, since y~a~ x dy = — |dy~a, we obtain 

(A) = - ({l- e - sP ^)y-^{~) r -/ sP m ge- sP ^y-^(~)dy)f(g)dg 
We have assumed a > 2, so lim^ofl — e~ sPmgy )y~a = 0. Letting i = sP m gy, we obtain 

i /-oo psP m gr- a 

(A) = i jf (-(l-e-^^)r 2 + ( S P m ^iy o e-V-^-My)/^)^ 

= -- + -^(sP m r-«) + -( s p m )^ / ^(r(i--)-r(i--, s p m ^- a ))/(^. 

22 2 Jq V Q: ct / 
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From the properties of Gamma function, we have 



r 











H 




- 2 r( 








a ' 


v a J 



r ( 1 - sP^r"^ = ( S P^r- Q )- 2 / a e- sP '^" a - -r ( --, sP m gr- 
\ a J a \ a 

so that (A) can be simplified as 

(A) = -^ + ^P m )^M(n--,sP m gr- a )-T(--))). 
2 a V V a a / / 

The evaluation of (5) is almost the same, and we have 

(5) = --(sWP f ^T(--)E 9 (g^. 

a a 

Substituting (A) and (B) into £/ o+ /„(s) and rearranging terms in Z m (T) hence leads to ©. 

Now we evaluate the mean achievable rate. Since for a positive random variable X, E[X] = 
J P{X > t}dt, we have 



P[ln(l + SINR)] 



,-?rA m r 2 



In 1 + 



OO /"00 

„-7rA m r 2 



I c + I f + ° 2 
Pmhr~ 



27rAr„rdr 



> t >dt2ir\ m rdr 



OO /"OO 

„-7rA m r 2 



OO /"OO 



OO /"OO 



JO 



™ f , / r m Air 
P{ ln(l + — jp- 
L i c + if 

/ F\h > —(I c +I f + (7 2 )(e* - 1) [dt27rA m rdr 

70 Pyn ' 

E(e^^ (Jc+ ^ +ff2)(e '- 1) )dt2 7 rA m rdr 
^(e*-i), /^r a (e*-l)- 



dt27rA rr ,rdr. 



(47) 



From the derivation of Z m {T) we have already obtained the expression of £>i c +i f 
When substituting into (l47l) . we arrive at (fTOl) and thus establish Theorem [IJ 



Appendix B 



Let r be the distance between a femtocell UE and its serving FAR Because the femtocell UEs 
are uniformly distributed in the circular coverage area of radius Rf of each FAP, the pdf of r is 
given by /(r) = 2r/R'j. 
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Denote by I' c and I'j the interference strengths from the macrocell tier and the femtocell tier 
respectively. Similar to the proof of Theorem [Q we have 

f R f 2r »Tr a <, 2 nTr a 
Z f (T) = l- *e- p f L^fLJL)^. (48) 

Jo ttf ' Pf 



Following the same procedure as in the derivation of (1451) . we get the Laplace transform of 

I' c + I'j as 

£ / , +/} ( S ) = ^ +/} (e-^ +/ /)) 

, poo poo 

= exp ( - 2it\' m J (1 - L g (sv- a WP m ))vdv - 2n\' f J (1 - L g (sw- a W 2 P f ))wdw 

= exp ( - vrsirXl - ^){X' m (WP m )^ + \' f (W 2 P f )%)E g (gty. (49) 

In the above, it is noteworthy that the interference from an FAP penetrates two walls thus the 
loss becoming W 2 instead of W. Substituting (|49| ) into (l48l) with v = r 2 , we get the SINR 
distribution in Theorem [2] Similar to (|47l) . the mean achievable rate is given by 

R f roc B r^ {et _ 1) / M r a (e*-1)\ 2r 



T ' = JJ e " £ '^l^V^«f dn <50) 

Substituting (|49|) into (l50l) . we get the mean achievable rate in Theorem |2] 
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TABLE I 
SYSTEM PARAMETERS 



Symbol 


Description 


Value 




PPPs defining the MBSs and FAPs 


N/A 


Am 


Density of MBSs 


0.00001 MBS/m 2 


A/ 


Density of FAPs 


0.0001 FAP/m 2 


Rf 


Radius of femtocell 


10m 


Aout 


Densities of outside nonsubscribers 


0.0001 user/m 2 


A s , Ai n 


Densities of subscribers and inside nonsubscribers 


0.015 user/m 2 


Pm, Pf 


Transmit power at MBS and FAP 


39dBm,13dBm 


M 


Number of subchannels at each access point 


20 


M a , M r 


Number of subchannels shared and reserved by each femtocell 


Not fixed 


a 


Path loss exponent 


4 


pm pj 
-'busy? -'busy 


Probabilities that a given subchannel is used by an MBS and an FAP 


Not fixed 


a' a'. 


PPPs defining the MBSs and FAPs that interfere a given subchannel 


N/A 


A m , A j 


Densities of MBSs and FAPs that interfere a given subchannel 


Not fixed 




Numbers of nonsubscribers that access a given MBS and a given FAP 


Not fixed 


U s 


Numbers of subscribers that access a given FAP 


Not fixed 


Tf,T m 


Mean achievable rates of femtocell UEs and macrocell UEs 


Not fixed 


7~S 5 ' 71 


Mean achievable rates of subscribers and nonsubscribers 


Not fixed 




Fig. 1. Downlink two-tier network model for hybrid access femtocell. 
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Fig. 2. The Voronoi network model. Each Voronoi cell is the coverage of a macrocell and each small circle stands for a 
femtocell. 



subchannels 



FAP 1 



FAP 2 



FAP 3 



Subchannels allocated to the nonsubscribers 
Subchannels reserved for the subscribers 



Fig. 3. Spectrum allocation in each hybrid access femtocell. All the M s shared subchannels are randomly selected in each 
FAP. 
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SINR Distribution 




SINR (dB) 



Fig. 4. Cdf of SINR for macrocell UEs and femtocell UEs when M s — 10. Other parameters are set as default. 




(a) Change the density of the FAPs A/. (b) Change the density of outside nonsubscribers A ou t. 

Fig. 5. Performance contrast between macrocell UEs and femtocell UEs with various densities of FAPs and outside 
nonsubscribers. 
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— © — Simulative Macrocell UEs 
■ © - Theoretical Macrocell UEs 
— H — Simulative Femtocell UEs 
Theoretical Femtocell UEs 



10 12 14 16 18 20 22 24 
The transmission power of the FAPs (dBm) 



5 5 - 



- Simulative Macrocell UEs 

- Theoretical Macrocell UEs 

- Simulative Femtocell UEs 
Theoretical Femtocell UEs 



30 35 40 45 50 

The transmission power of the MBSs (dbm) 



(a) Change P f while P m = 39dBm. (b) Change P m while Pj = 13dBm. 

Fig. 6. Performance contrast between macrocell UEs and femtocell UEs with various transmission powers. 
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■9 — Simulative Nonsubscribers 
© - Theoretical Nonsubscribers 
■a — Simulative Subscribers 
i - Theoretical Subscribers 



2 4 6 8 10 12 14 16 18 20 
The number of shared subchannels in each FAP M 



Fig. 7. Performance of nonsubscribers and subscribers with various numbers of shared subchannals M s in an FAP. Other 
parameters are set as default. 
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Simulative Nonsubscribers 
Theoretical Nonsubscribers 
Simulative Subscribers 
Theoretical Subscribers 



2 4 - 




2 4 6 8 10 12 14 16 18 20 
The number of shared subchannels in each FAP M 



- Simulative Nonsubscribers 

- Theoretical Nonsubscribers 

- Simulative Subscribers 
Theoretical Subscribers 



2 4 6 8 10 12 14 16 18 20 
The number ot shared subchannels in each FAP M 



(a) A s = 0.003 and A in = 0.027. 



(b) X a = 0.027 and A in = 0.003. 



Fig. 8. Performance contrast between nonsubscribers and subscribers with different proportions of inside nonsubscribers and 
subscribers. 
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